(Leave  blank) 


AD 


Award  Number: 

W81XWH-07-1-0149 

TITLE : 

Exploiting  a  Molecular  Gleason  Grade  for  Prostate  Cancer  Therapy 


PRINCIPAL  INVESTIGATOR: 
Peter  S.  Nelson,  M.  D. 


CONTRACTING  ORGANIZATION: 

Fred  Hutchinson  Cancer  Research  Center 
Seattle,  WA  98109 

REPORT  DATE: 

March  2009 


TYPE  OF  REPORT: 
Annual 


PREPARED  FOR:  U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


DISTRIBUTION  STATEMENT: 

X  Approved  for  public  release;  distribution  unlimited 

Distribution  limited  to  U.S.  Government  agencies  only; 
report  contains  proprietary  information 


The  views,  opinions  and/or  findings  contained  in  this  report  are 
those  of  the  author  (s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy  or  decision  unless  so 
designated  by  other  documentation. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

03/14/2009  Annual 

3.  DATES  COVERED  (From  -  To) 

02/15/2008-02/14/2  0  09 

4.  TITLE  AND  SUBTITLE 

Exploiting  a  Molecular  Gleason  Grade  for  Prostate  Cancer  Therapy 

5a.  CONTRACT  NUMBER 

W81XWH-07-1-014  9 

5b.  GRANT  NUMBER 

PC0  605  95 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Peter  S.  Nelson,  M.D. 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Fred  Hutchinson  Cancer  Research  Center 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

1100  Fairview  Ave  N 
PO  Box  19024 
J6-500 

Seattle,  WA  98109-1024 

9.  SPONSORING  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Medical  Research 
And  Materiel  Command 
Fort  Detrick,  MD  21702-5012 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  Purpose  of  this  proposal  is  to  exploit  a  molecular  correlate  of  the  Gleason  grading 
system  for  prostate  carcinoma  in  order  to:  a)  develop  improved  outcome  predictors;  and  b) 
identify  therapeutic  strategies.  During  this  project  period  we  have  evaluated  in  excess  of  25 
prostate  cancer  antigens  by  Western  blot  and  tissue-based  assays  and  identified  3  with 
detectable  levels  in  the  plasma.  We  also  determined  that  developmental  genes  involved  with 
invasive  cellular  behavior  exhibit  associations  with  clinical  outcomes.  While  these  results 
are  important,  the  further  exploitation  of  these  findings  is  limited  by  the  lack  of  specific 
antibodies  suitable  for  tissue  (IHC)  and  ELISA-based  studies  of  protein  localization  and 
quantitation.  Ongoing  studies  are  focused  on  developing  alternative  strategies  for  protein 
quantitation  in  order  to  complete  the  proposal  objectives. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

USAMRMC 

a.  REPORT 

U 

b.  ABSTRACT 

U 

c.  THIS  PAGE 

U 

UU 

18 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


Table  of  Contents 


Page 

Introduction .  4 

Body .  4-8 

Key  Research  Accomplishments .  8 

Reportable  Outcomes .  9 

Conclusion .  9 

References .  9 


Appendices 


9 


Nelson-PC060595-Progress  Report  02/09-Year  2 


INTRODUCTION 

This  proposal  is  designed  to  exploit  a  molecular  correlate  of  the  Gleason  grading  system  for 
prostate  carcinoma  in  order  to:  a)  develop  improved  outcome  predictors;  and  b)  identify  thera¬ 
peutic  strategies  targeted  toward  features  unique  to  aggressive  cancers.  We  hypothesized  that  the 
specific  molecular  features  that  underlie  prostate  cancer  grades  define  the  capacity  for  tumor  cell 
invasion  and  dissemination  (progression)  and  may  represent  unique  diagnostic  markers  and  tar¬ 
gets  for  therapeutic  intervention. 

The  aims  of  the  proposal  remain  unchanged:  Aim  #1:  To  compare  the  power  of  molecular  versus 
histological  Gleason  categories  for  outcome  predictions  in  the  context  of  PSA  relapse  and  pros¬ 
tate  cancer-specific  mortality.  Aim  #2:  Determine  if  grade-associated  differences  in  prostate  can¬ 
cer  protein  expression  are  reflected  by  levels  of  cognate  serum  proteins.  #3  Write  final  report. 
(Note  the  original  Aim  3  involving  animal  studies  of  altering  prostate  cancer  grade-associated 
functions  was  deleted  due  to  recommendations  by  the  reviewers). 

Disease  relevance:  Through  comprehensive  studies  of  genome  and  gene  expression  alterations,  it 
is  clear  that  prostate  cancers  are  profoundly  heterogeneous,  both  at  the  molecular  and  clinical 
level.  In  this  context,  the  successful  completion  of  this  proposal  has  the  potential  to  greatly  im¬ 
prove  clinical  decision  making  by  improving  the  accuracy  of  predicting  which  cancers  may  best 
be  ‘treated’  by  observation  rather  than  primary  therapy,  or  defining  those  cancers  that  should 
have  additional  systemic  therapy— in  addition  to  local  therapy— by  virtue  of  their  high  malignant 
potential.  Also  in  the  context  of  clinical  care,  several  of  the  key  nodes  distinguishing  low  grade 
from  high-grade  cancer  feature  metabolic  components  that  currently  can  be  targeted  by  FDA- 
approved  drugs  originally  designed  for  the  treatment  of  diseases  other  than  neoplasia  (e.g. 
MAO  A  inhibitors). 

BODY 

The  following  summarizes  the  technical  objectives  for  the  proposal  and  the  work  accomplished 
during  the  12-month  interval  between  the  Year  1  Progress  Report  (02/08)  and  the  present  Year  2 
Progress  Report  (02/09).  Work  done  during  Year  1  (02/07-02/08)  and  reported  last  year  is  in¬ 
cluded  in  parentheses  (..)  as  a  frame  of  reference. 

D.l. Technical  objective  1:  To  compare  the  predictive  power  of  molecular  versus  histological 
Gleason  categories  for  outcome  predictions  in  the  context  of  PSA  relapse  and  prostate  cancer- 
specific  mortality.  (Months  1-24). 

Objective  la.  Antibody  acquisition  and  evaluation. 

Task  1 :  purchase  antibodies  recognizing  grade-determinant  proteins  (months  1-12). 

(Year  1:  To  date,  we  have  purchased  (or  acquired)  antibodies  recognizing;  TMPRSS2,  MAOA, 
DAD1,  ERG,  Jagged,  p63,  AMACR,  MUC1,  FLNA,  ALSCR2,  CCNG2,  FLH2,  GSTMU1, 
PC4,  RSK2,  and  SMS.  The  majority  of  these  have  previously  been  used  in  immunohisto- 
chemical  studies,  though  not  in  prostate  cancer.) 

Year  2:  We  have  purchased  an  additional  25  antibodies  for  the  evaluation  of  prostate  cancer  as¬ 
sociated  gene/protein  expression:  Postn,  Arfl,  Arf2,  Cite,  Famc2,  Msn,  Nrpl,  Ppplcc, 
Psma7,  Tpm3,  Tubala,  Tubalb,  Tuba4a,  AZGP1,  PTPRF,  ADAM7,  ADAM9,  AD  AMTS  5, 
AFC  AM,  CXCF12,  EFNA1,  HSPD1,  JAG2,  NOTCH3,  and  STOMF2. 

Task  2:  evaluate  each  antibody  using  a  semi-quantitative  immunoblot  (months  1-12). 
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/Year  1:  We  performed  semi-quantitative  immunoblots  for  approximately  1/2  of  these  antibod¬ 
ies  with  the  specific  intent  of  demonstrating  a  specific  band.  However,  we  have  found  that 
this  approach  is  not  definitive,  and  does  not  provide  additional  data  in  our  hands  due  to  the 
profound  heterogeneity  in  the  tissue,  and  highly  variable  expression  in  cell  lines.) 

Year  2:  We  performed  semi-quantitative  immunoblots  for  25  additional  antibodies  (see  Task  1 
above  for  Year  2)  recognizing  that  despite  the  limitations  of  the  approach,  when  moving  for¬ 
ward  to  tissue  staining,  it  is  important  to  have  an  assessment  of  antibody  specificity  when  in¬ 
terpreting  immunostains— e.g.  avoid  a  false-positive  signal.  While  very  few  of  these  com¬ 
mercial  antibodies  provide  a  single  band  of  expected  protein  size,  several  resulted  in  quite 
promising  distinction  between  benign  and  cancerous  prostate  {Figure  1,  Below  Left,  PTPRF 
shows  elevated  expression  in  prostate  and  breast  carcinoma,  Figure  2,  Right,  JAG2  shows 
elevated  expression  in  prostate  cancer  and  in  several  prostate  cancer  cell  lines,  and  condi¬ 
tioned  medium  from  these  lines,  indicating  protein  secretion),  and  several  also  detected  an 
appropriate  protein  in  serum  from  tumor  bearing  xenografts  (see  Figure  3,  below). 


PTPRF  (Protein  Tyrosine  Phosphatase.  Receptor-type  F) 
Western  Blot 

(Predicted  molecular  weight:  Cleaved  extracellular  domain  139  kOa) 


1  2345678 


Lanes  (30  M^lane): 

1.  Human  universal 

2.  Human  Liver 

X  Human  Prostate:  Normal 

4.  Human  Prostate:  Tumor 

5.  Human  Breast:  Normal 
&  Human  Breast:  Tumor 
7.  Mouse  Universal 

&  Mouse  Liver 


1  f  Goat  and-PTPRR  i  :500  ( HAD  Systems) 

2f  HRP-Conjugaied  Mouse  ant f  Goat  1:10.000  ( Pierce.' Thermo) 

2  mm  exposure  (TW-06110G) 


JAG2  (jagged  2 ) 

Western  Blot 

(Predicted  molecular  weight:  129.2  kDa;  -13X4  kOa) 


1  *:  RaDOit  amt-JAG2  i  :i  .000  (Cell  Signaling  Technology) 

2*  HRP-Coni'jqated  Mouse  anH-RaOblt.  l:iaooo (Pierce Thermo) 
3  min  exposure  (JT-OS  UOffl 


Figure  1.  Elevated  expression  of  PTPRF  protein  in 
prostate  cancer  tissue. 

Figure  2.  Elevated  expression  of  JAG2  protein  in 
prostate  cancer  tissue. 

Task  3:  optimize  IHC  staining  using  a  fixation  TMA,  antibody  dilutions  and  antibody  re¬ 
trieval  methods  (months  1-12). 

/Year  1:  We  have  optimized  titering  and  conditions  for  the  following  antibodies  against  a 
small  panel  of  benign  and  neoplastic  prostate  cases:  TMPRSS2,  MAOA,  DAD1,  ERG,  Jagged, 
p63,  AMACR,  MUC1,  FLNA,  ALSCR2,  CCNG2,  FLH2,  GSTMU1,  PC4,  RSK2,  and  SMS— 
see  reportable  outcomes,  Datta  et  al  (2007). 

Year  2:  We  have  continued  optimization  of  antibody  immunostaining  for  additional  antibod¬ 
ies  associating  with  grade/prostate  cancer  outcome.  These  have  included  ABP280L,  ITGAS, 
PSA,  MIB1,  Mucl,  Jagged  1,  AMACR,  CDK7,  MTA1  and  SKC4A1AP. 

Objective  lb.  IHC  analysis/confirmation  of  protein  expression  and  Gleason  pattern 

Task  4:  compare  protein  expression  patterns  relative  to  transcript  measures  by  microar¬ 
ray  (months  12-16). 

(Year  1:  We  have  initiated  these  studies.  We  have  re-organized  grade-defining  genes  present 
on  our  original  microarray  studies  into  a  format  that  can  be  directly  compared  to  protein  expres¬ 
sion  levels  determined  by  IHC.  We  are  now  in  the  process  of  quantitating  the  protein  expression 
across  Gleason-appropriate  prostate  cancer  cases,  and  will  import  the  data  directly  for  compari- 
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son.  We  have  completed  a  thorough  analysis  of  the  expression  of  TMPRSS2,  a  grade-associated 
protein.  Figure  1  demonstrates  TMPRSS2  expression  in  benign  and  cancer  epithelium  and  shows 
mislocalization  of  TMPRSS2  in  cancer  cells.  Figure  2  demonstrates  Gleason  Grade-associated 
TMPRSS2  expression.) 

Year  2:  We  are  continuing  these  studies.  To  date,  very  few  of  the  antibodies  have  yielded 
consistent  transcript  and  protein  correlations.  TMPRSS2  is  an  exception  where  we  found  clear 
correlations  between  transcript  and  protein  upregulation  in  prostate  cancer  that  also  associated 
with  Gleason  grade  (see  Reportable  Outcomes  from  Year  1,  Lucas  et  al).  A  component  of  this 
discordance  is  likely  due  to  the  protein  localization  in  different  cellular  compartments  (e.g. 
membrane  vs  nucleus).  Another  component  of  this  discordance  may  lie  in  the  lack  of  antibody 
specificity.  A  third  component  of  this  discordance  may  simply  reflect  lack  of  clear  concordance 
between  transcript  and  protein  levels  in  cells/tissues. 

Objective  lc.  1HC  analysis  of  prostate  cancer  cohorts  with  outcomes  reflected  by  PSA  relapse 
and  mortality. 

Task  5:  Stain  and  read  Gelman  TMAs:  20  antibodies  (months  16-22) 

Year  2:  We  completed  ‘outcomes’  studies  of  the  TMPRSS2  antibody  on  the  Gelman  TMA 
and  found  no  correlation  with  relapse  that  was  independent  of  Gleason  Grade  (see  Reportable 
Outcomes  for  Year  1,  Lucas  et  al).  In  the  process  of  conducting  these  studies,  we  determined  that 
this  TMA  is  not  optimal  for  completing  all  of  the  planned  studies  due  to  the  limited  numbers  of 
these  arrays  that  will  actually  be  available.  Thus,  we  have  initiated  the  construction  additional 
TMAs  with  outcomes.  We  have  developed  a  network  of  collaborating  institutions  (5  different 
sites)  that  will  each  construct  outcomes  TMAs  suitable  for  cross-comparison  of  IHC  results  in 
the  context  of  outcomes  associations.  We  have  identified  200  cases  within  our  Center  that  will  be 
used  and  construction  of  this  TMA  has  commenced  with  planned  completion  by  6/09.  The 
planned  IHC  studies  for  Task  5  will  commence  at  that  time.  We  have  constructed  a  50  case  out¬ 
comes  TMA  and  have  used  this  resource  to  evaluate  outcomes  of  4  markers  to  date,  (see  Task  7). 

Task  6:  Stain  and  read  Stanford  TMAs:  20  antibodies  (months  16-22).  The  construction 
of  the  Stanford  TMA  was  delayed  due  to  logistics  of  acquiring  all  sample  blocks  by  Dr.  Stan¬ 
ford’s  group.  These  blocks  have  now  been  collected  and  the  construction  of  this  TMA  has  com¬ 
menced  with  an  anticipated  completion  date  of  07/09. 

Year  2: 

Task  7:  Determine  statistical  associations  with  outcomes  (months  22-24).  We  found  no  in¬ 
dependent  correlations  with  TMPRSS2  expression  and  clinical  outcomes  using  the  original  Gel¬ 
man  TMA.  We  have  completed  a  study  of  4  grade/outcome-associated  antibodies  (ITGAS,  PSA, 
MIB1,  MTA1)  and  found  no  correlations  with  outcomes  using  the  50-case  outcomes  TMA. 
These  will  be  repeated  using  the  200  case  outcomes  TMAs  and  the  Stanford  TMA. 

In  a  companion  study  that  involved  assessing  molecular  correlates  of  prostate  development, 
we  identified  genes  associated  with  the  branching  morphogenesis/invasive  component  of  devel¬ 
opment  and  features  of  prostate  carcinoma.  Several  of  these  genes  were  Gleason  grade- 
associated  genes.  We  examined  existing  prostate  cancer  datasets  with  annotated  clinical  out¬ 
comes,  and  determined  that  several  of  these  developmental  genes  associated  with  prostate  cancer 
relapse  following  primary  therapy  (see  Figure  4,  in  the  appended  manuscript  provided  in  Report- 
able  Outcomes,  Pritchard  et  al).  This  study  provides  additional  candidates  to  evaluate  in  the  con¬ 
text  of  grade/outcome  markers  for  prostate  cancer  that  is  the  focus  of  the  present  application,  and 
we  are  currently  obtaining  antibodies  for  verification  studies  based  on  tissue  protein  expression. 

Task  8:  Refine  antibody/protein  list  to  minimal  redundant  set.  (months  22-24) 

Year  2:  This  task  is  pending  the  completion  of  the  TMA  studies  above  in  Tasks  6-7. 
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D.2.  Technical  ob  jective  2:  Determine  if  grade-associated  differences  in  prostate  cancer  pro¬ 
tein  expression  are  reflected  by  levels  of  serum  proteins  (months  3-36). 

Objective  2a.  Western  Analysis  for  Antibody  Q/C. 

Task  9:  prepare  Western  blots  of  serum  proteins  (months  3-6). 

(Year  1:  We  have  acquired  a  panel  of  (anonymized)  human  serum  protein  samples  that  span 
a  spectrum  of  a)  absence  of  prostate  cancer — biopsy  proven;  b)  low  grade  prostate  cancer;  c) 
high  grade  prostate  cancer;  d)  metastatic  prostate  cancer.  The  quality  of  the  samples  has  been 
verified  using  Western  analysis  for  abundant  and  low  abundant  proteins. 

Task  10:  determine  specificity  of  immunoreactivity  and  semi-quantitation  (months  7-12). 

(Year  1:  We  have  completed  Western  analysis  (blots)  for  12  proteins/antibodies.  Of  these, 
four  antibodies  produced  patterns  indicating  poor  specificity,  with  multiple  bands  present.  Of  the 
remaining  8,  three  did  not  provide  satisfactory  detection,  and  these  are  being  re-evaluated,  or  an¬ 
other  antibody  source  is  being  evaluated.  The  remaining  five  demonstrating  good  specificity  and 
we  are  scaling  up  a  larger  sample  set  to  assess  overall  cancer/benign/grade  evaluation.  One  pro¬ 
tein,  osteopontin,  is  capable  of  distinguishing  metastatic  cancer,  but  not  early  stage  cancer,  rela¬ 
tive  to  individuals  without  cancer.) 

Year  2:  We  have  completed  the  serum-based  analysis  of  an  additional  25  antibodies  recog¬ 
nizing  prostate  cancer  antigens  (see 
Objective  la,  Task  1,  Year  2).  Of  these,  10 
recognized  an  appropriate  band  in 
conditioned  medium  from  prostate  cells 
and  3  recognized  an  appropriate  band  in 
the  plasma  of  mice  bearing  prostate 
tumors  (see  Figure  3).  Note,  we  have 
opted  to  use  a  stepwise  approach  for 
evaluating  tumor-associated  plasma 
proteins  to  save  valuable  human  plasma. 

(The  murine  plasma  samples  were 
available  from  other  ongoing  studies  and 
thus  were  not  a  cost  to  this  project).  The 
next  phase  will  involve  the  direct 
evaluation  of  these  candidates  in  human 
samples  by  Western  blot  and  ELISA. 

Objective  2b.  ELISA 

Task  11:  Prepare  assay  plates  (months  6-10). 

(Year  1:  We  have  prepared  ELISA  plates  for  three  antibodies  combinations  to  date.  Data  for 
our  studies  of  osteopontin  are  shown  in  Figure  1.  We  have  established  a  system  to  avoid  non¬ 
specific  detection  for  these  proteins.  For  the  remainder  of  candidates,  we  are  attempting  to  iden¬ 
tify  a  second  antibody  as  most  ELISAs  rely  on  two  antibodies  for  specificity.  If  we  are  unable  to 
identify  a  second  antibody,  we  will  proceed  with  a  capture-ELISA  or  move  toward  other  novel 
approaches  such  as  SISCAPA  (mass  spectrometry-based  assay — see  task  16)). 

Year  2:  We  have  been  unable  to  identify  suitable  secondary  antibodies  for  the  development 
of  sandwich  ELISAs  for  the  candidate  protein  biomarkers  (with  the  exception  of  osteopontin). 
The  capture  ELISAs  have  not  been  of  sufficient  sensitivity  or  specificity  to  detect  protein  levels. 
We  are  continuing  to  evaluate  new  ELISAs  developed  by  commercial  manufacturers  to  our  can- 


PSMA7  (proteasome  (prosome.  macropain)  subunit  alpha  type,  7) 
Western  Blot 

(Predicted  motecuiar  weight:  -.27.9  kOa) 


191 - ■ 

97  — ■ 

64 - , 

51 - T 

30 - ■ 

20  —  .o 


123456789 


UZ3  [==] 


Lanes <20  up  lane): 

1.  22Rvi  control  Pool  i :  Depleted  Plasma' 
2  22Rvl  Tumor  Pool  1 1  Depieed  Plasma* 

1  22Rvi  Tumor  Pool  2:  Depleted  Plasma* 

4.  LNCaP  Control  Pool  1 :  Depleted  Plasma* 

6.  LNCaP  Tumor  Pool  1 :  Depleted  Plasma* 
&  LNCaP  Tumor  Pool  2  Depleed  Plasma* 

7.  Blank  Weil 

&  22Rvi  xenograft  Tissue  Lysate  Pool 
fl.  LNCaP  Xenograft  Tissue  Lysate  Pool 

•pools  of  3  animals  each 


1  RaPOlt  antf-20S  Proteasome  07:o&  1 :200  (Santa  Cruz) 

2^  HRP  Conjugated  Mouse  antl-RabNt  1:10.000 (Pierce,1 THermo) 
3  man  exposure  (TW) 


Figure  3.  Expression  of  PSMA7  in  the  plasma  of  mice  bear¬ 
ing  prostate  cancer  xenografts  (boxes).  Note,  absence  of 
PSMA7  in  serum  of  control  mice  (lanes  1  and  4). 
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didate  proteins,  and  we  have  initiated  the  development  of  new  antibodies  with  a  collaborator,  Dr. 
Brad  Nelson,  at  the  University  of  British  Columbia.  We  will  evaluate  these  antibodies  for  sensi¬ 
tivity/specificity  as  they  become  available  early  in  Year  3. 

Task  12:  Run  q/c  with  recombinant  protein  standards,  (months  10-11). 

(Year  1:  For  two  proteins,  we  have  run  q/c  for  two  proteins  and  these  have  passed  our  met¬ 
rics.) 

Year  2:  This  has  proven  to  be  a  significant  challenge.  We  have  run  q/c  for  additional  pro¬ 
teins  and  no  ELISAs  have  passed  q/c  metrics.  For  this  reason,  we  have  opted  to  develop  new  an¬ 
tibodies  as  described  above  for  Task  11,  and  we  have  moved  forward  with  SISCAPA  analyses 
(see  Task  16).  We  have  also  initiated  studies  using  single-antibody  assays  by  the  Luminex  (bead- 
based)  system. 

Task  13:  Analyze  control  and  disease  serum  samples  by  ELISA  (months  10-14). 

(Year  1:  We  have  recently  run  serum  samples  representing  low  grade  and  high-grade  cancer 
as  well  as  benign  and  metastasis.  Data  analysis  is  in  progress.) 

Year  2:  We  completed  the  sample  analysis  for  Osteopontin  and  a  manuscript  detailing  these 
results  is  in  preparation  (note,  details  of  these  results  were  shown  in  the  Year  2  progress  report). 
Task  14:  Repeat  Tasks  10-13  for  two  additional  proteins  (months  14-24)-pending 
Year  2:  Due  to  difficulties  obtaining  two  high-quality  antibodies/candidate  protein  (required 
for  capture  ELISAs),  we  have  yet  to  develop  ELISAs  for  the  other  promising  protein  candidates 
(e.g.  PSMA7-Figure  3).  We  are  currently  evaluating  the  Luminex  system  for  the  utilization  of 
single-antibodies  and  continue  to  assess  commercial  inventories  of  new  ELISAs  being  marketed. 
We  have  also  begun  studies  using  SISCAPA  (see  Objective  2c,  below). 

Task  15:  Repeat  Tasks  10-13  for  two  additional  proteins  (months  24-34)-pending 

Objective  2c.  SISCAPA 

Task  16:  Explore  the  utility  of  SISCAPA  (Stable  Isotope  Capture  by  Anti-Peptide  Antibodies) 
as  an  alternative  (and  improvement)  to  ELISA-based  assays,  (months  30-35). 

Year  2:  Due  to  the  challenges  we  have  observed  in  obtaining  multiple  high-quality  (specific) 
antibody  reagents,  we  have  started  evaluating  the  SISCAPA  approach  for  protein  quantitation. 
This  involves  the  development  of  mass  spectrometry  peptide  sprectra  for  the  candidate  proteins. 
We  have  now  identified  highly  specific  spectra  for  peptides  corresponding  to  16  of  the  25  candi¬ 
date  proteins  described  in  Specific  Aim  la,  Year  2.  We  are  currently  determining  if  these  pep¬ 
tides  can  be  observed  by  mass  spectrometry  in  the  serum  of  tumor-bearing  hosts. 

D.3.  Technical  objective  3:  Final  Report:  Complete  data  analyses,  compile  accomplishments 
and  reportable  outcomes  and  write  final  project  report  (Months  35-36)-pending. 

KEY  RESEARCH  ACCOMPLISHMENTS  (Year  2) 

•  We  completed  the  acquisition  and  preliminary  analysis  of  25  additional  antibodies  that  rec¬ 
ognize  a  subset  of  Gleason  Grade-associated  molecular  changes  found  in  human  prostate 
cancer.  Several  were  capable  of  distinguishing  benign  from  malignant  prostate  and  3  recog¬ 
nized  a  protein  of  appropriate  size  in  the  plasma  of  tumor-bearing  mice.  The  development  of 
quantitative  detection  methods  for  these  proteins  is  in  progress  (e.g.  ELISA,  SISCAPA). 

•  We  integrated  a  study  of  prostate  development  and  prostate  carcinoma  to  identify  a  cohort  of 
prostate  cancer  genes  that  associate  with  clinical  outcomes,  a  subset  of  which  include  Glea¬ 
son-associated  genes  (see  Reportable  Outcomes).  Characterization  of  protein-based  quantita¬ 
tion  is  in  progress. 
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•  We  assembled  the  requisite  tissue  resources  to  construct  tissue  microarrays  (TMAs)  of  suit¬ 
able  quality  and  quantity  (number  of  tissue  sections  to  allow  for  the  planned  antibody  assays) 
to  enable  outcome  studies  of  candidate  Gleason  grade-associated  proteins. 

REPORTABLE  OUTCOMES 

Year  1: 

Datta  MW,  True  LD,  Nelson  PS,  Amin  MB.  The  role  of  tissue  microarrays  in  prostate  cancer 
biomarker  discovery.  (2007)  Adv  Anat  Pr?//ro/.Nov;14(6):408-18. 

Lucas  J,  True  L,  Hawley  S,  Matsumura  M,  Morrissey  C,  Vessella  R,  Nelson  P.  (2008)  The  an¬ 
drogen-regulated  type  II  serine  protease  TMPRSS2  is  differentially  expressed  and  mislocal- 
ized  in  prostate  adenocarcinoma.  J  Pathol.  215:118-125. 

Year  2: 

Pritchard  C,  Mecham  B,  Dumpit  R,  Coleman  I,  Bhattacharjee  M,  Chen  Q,  Sikes  RA,  Nelson  PS. 
(2009).  Conserved  gene  expression  programs  integrate  mammalian  prostate  development  and 
tumorigenesis.  Cancer  Res.  2009  Mar  l;69(5):1739-47. 

CONCLUSIONS 

We  have  demonstrated  that  cancer  and  grade-associated  mRNA  abundance  levels  are  associated 
with  corresponding  protein  alterations  for  a  subset  (but  not  all)  of  candidates.  The  major  limita¬ 
tion  to  our  studies  at  this  juncture  is  the  lack  of  availability  of  high  quality  antibody  reagents  that 
are  essential  for  developing  a  protein-based  (IHC  and  serum)  correlate  to  our  transcript  studies  of 
Gleason  grade.  This  is  not  unique  to  our  studies,  as  high-quality  antibody  reagents  are  a  limiting 
factor  for  many  studies  of  protein  correlates  of  gene  expression  in  cells  and  tissues.  However,  we 
have  several  promising  candidates  based  on  our  studies  accomplished  to  date  with  available  anti¬ 
bodies  and  have  embarked  on  several  alternative  strategies  that  we  anticipate  will  allow  for  the 
completion  of  the  project.  These  include:  (i)  development  of  bead-based  protein  assays  that  allow 
for  the  use  of  a  single,  rather  than  2,  antibodies;  (ii)  development  of  SISCAPA-based  approaches 
for  mass  spectrometry-based  protein  quantitation;  (iii)  development  of  new  antibodies. 
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Abstract 

Studies  centered  at  the  intersection  of  embryogenesis  and 
carcinogenesis  have  identified  striking  parallels  involving 
signaling  pathways  that  modulate  both  developmental  and 
neoplastic  processes.  In  the  prostate,  reciprocal  interactions 
between  epithelium  and  stroma  are  known  to  influence 
neoplasia  and  also  exert  morphogenic  effects  via  the 
urogenital  sinus  mesenchyme.  In  this  study,  we  sought  to 
determine  molecular  relationships  between  aspects  of  normal 
prostate  development  and  prostate  carcinogenesis.  We  first 
characterized  the  gene  expression  program  associated  with 
key  points  of  murine  prostate  organogenesis  spanning  the 
initial  in  utero  induction  of  prostate  budding  through 
maturity.  We  identified  a  highly  reproducible  temporal 
program  of  gene  expression  that  partitioned  according  to 
the  broad  developmental  stages  of  prostate  induction, 
branching  morphogenesis,  and  secretory  differentiation. 
Comparisons  of  gene  expression  profiles  of  murine  prostate 
cancers  arising  in  the  context  of  genetically  engineered 
alterations  in  the  Pten  tumor  suppressor  and  Myc  oncogene 
identified  significant  associations  between  the  profile  of 
branching  morphogenesis  and  both  cancer  models.  Further, 
the  expression  of  genes  comprising  the  branching  morpho¬ 
genesis  program,  such  as  PRDX4,  SLC43A1,  and  DNMT3A,  was 
significantly  altered  in  human  neoplastic  prostate  epithelium. 
These  results  indicate  that  components  of  normal  develop¬ 
mental  processes  are  active  in  prostate  neoplasia  and  provide 
further  rationale  for  exploiting  molecular  features  of  organ¬ 
ogenesis  to  understand  cancer  phenotypes.  [Cancer  Res 
2009;69(5):  1739-47 ] 

Introduction 

Studies  involving  normal  developmental  processes  have  revealed 
important  parallels  with  carcinogenesis  that  involve  key  signaling 
mechanisms  controlling  the  three-dimensional  growth  and  orga¬ 
nization  of  tissues  (1,  2).  Organogenesis  is  a  complex  process 
involving  proliferation,  pattern  specification,  and  cellular  differen¬ 
tiation  orchestrated  by  an  evolving  transcriptional  program  (3). 
Many  highly  conserved  pathways  instrumental  in  dictating  ordered 
organ  and  organismal  morphogenesis,  originally  defined  in  model 
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organisms  such  as  Drosophila  melanogaster  and  Caenorhabditis 
elegans,  have  been  found  to  be  altered  in  human  cancers.  Examples 
include  networks  involving  Wnt/adenomatous  polyposis  coli/ 
catenins,  Notch/Delta/Jagged,  fibroblast  growth  factor,  epidermal 
growth  factor,  transforming  growth  factor  [1  (TGF(3>)/Smad,  and 
Hedgehog/Patched/Smoothened.  Importantly,  information  trans¬ 
mitted  via  pathways  controlled  by  these  molecular  interactions 
dictates  cellular  behaviors  beyond  mitogenic  responses  to  include 
positional  sense,  differentiation,  invasion,  motility,  the  production 
of  matrix  components,  and  synthesis  of  autocrine  and  paracrine 
signaling  molecules. 

Key  features  of  normal  prostate  organogenesis  involve  character¬ 
istics  that  are  also  hallmarks  of  prostate  neoplasia,  including  a 
dependence  on  hormonal  signaling,  severing  of  cell-cell  contacts, 
invasion  of  epithelium  into  the  organ  microenvironment,  cell 
migration,  reestablishment  of  cell  contacts,  and  the  development  of 
new  blood  vessel  networks  (4).  The  prostate  gland  is  an 
endodermal  derivative  of  the  hindgut  first  formed  in  late  fetal  life 
when  androgen  produced  by  the  testis  induces  urogenital  sinus 
(UGS)  epithelial  invasion  into  the  mesodermally  derived  UGS 
mesenchyme  (5).  The  vast  majority  of  work  detailing  prostate 
developmental  processes  has  involved  rodents  in  which  the  first 
prostate  buds  are  visible  at  day  17  of  embryogenesis.  Branching 
morphogenesis  shortly  follows  this  inductive  phase  and  proceeds 
through  the  first  15  days  of  postnatal  life  (6).  Androgen  levels 
steeply  rise  at  puberty  (25-30  days  postnatal  in  the  mouse), 
resulting  in  prostate  growth  and  terminal  secretory  differentiation 
that  is  complete  by  ~  45  days  postnatal.  Thus,  the  major  events  of 
mouse  prostate  development  can  be  summarized  in  three  broad 
steps  that  comprise  prostate  induction,  branching  morphogenesis, 
and  secretory  differentiation. 

To  date,  the  fundamental  molecular  processes  mediating  the 
malignant  phenotypes  of  prostate  cancer  cells  remain  poorly 
defined.  Because  only  limited  temporal  information  can  be  gained 
from  studies  of  any  discrete  focus  of  malignancy,  we  reasoned  that 
systematically  evaluating  normal  cellular  processes  that  share 
features  with  prostate  carcinogenesis  may  provide  insights  into 
additional  networks,  pathways,  or  individual  molecular  interac¬ 
tions  that  contribute  to  neoplastic  growth.  Prostate  organogenesis 
and  carcinogenesis  both  exhibit  a  dependence  on  androgenic 
hormones,  and  each  is  influenced  by  a  complex  cross-talk  of 
paracrine  factors  operating  between  epithelium  and  stroma 
(mesenchyme;  ref.  5).  Further,  alterations  in  key  developmental 
signaling  nodes,  such  as  the  Sonic  Hedgehog  and  Notch  networks, 
exhibit  reproducible  alterations  in  prostate  carcinomas  (7-9). 
Based  on  these  findings,  we  sought  to  determine  the  relationships 
between  the  global  genetic  programs  associated  with  the  course  of 
prostate  organogenesis  and  those  found  to  be  influenced  by 
oncogenic  pathways  leading  to  invasive  cancer.  Herein,  we  detail 
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the  characterization  of  transcriptional  profiles  corresponding  to 
intervals  of  normal  prostate  organogenesis,  define  their  relation¬ 
ships  to  genetically  engineered  mouse  models  of  prostate 
carcinoma,  and  investigate  associations  with  human  prostate 
cancer  phenotypes. 

Materials  and  Methods 

Mouse  prostate  tissue  dissection,  RNA  isolation,  and  RNA  amplifi¬ 
cation.  All  mouse  studies  were  performed  in  accordance  with  Institutional 
Animal  Care  and  Use  Committee-approved  protocols.  Whole  male  UGS 
(E14.5,  E15.5,  E16.5,  E17.5,  and  E18.5)  or  separated  prostate  lobes  (P7,  P30, 
and  P90)  were  dissected  from  C57BL6/J  mice  and  snap  frozen  in  liquid 
nitrogen.  For  each  biological  replication,  we  pooled  3  to  10  mice 
representing  one  or  two  litters.  RNA  from  pools  of  UGS  or  specific  prostate 
lobes  (vp,  ap,  and  dip)  was  prepared  using  the  Qiagen  RNeasy  Mini  kit.  We 
included  an  on-column  DNasel  treatment  to  remove  contaminating  DNA. 
Before  RNA  amplification,  we  combined  equal  quantities  of  RNA  from  vp, 
ap,  and  dip  for  the  postnatal  prostate  samples.  We  amplified  1  pg  of  total 
RNA  from  each  sample  through  one  round  using  the  Arcturus  RiboAmp  kit. 
For  the  E14.5  UGS  reference  sample,  a  second  round  of  amplification  was 
done  to  provide  enough  RNA  for  all  microarrays.  Quantitative  reverse 
transcription-PCR  (RT-PCR)  showed  no  significant  difference  in  the  relative 
expression  of  two  genes  ( Plzf  and  Pdgfd )  between  unamplified  and 
amplified  RNA. 

Microarray  analyses.  cDNA  microarrays  enhanced  for  genes  expressed 
in  the  developing  mouse  prostate  were  prepared  as  previously  described 
(see  Supplementary  Data;  ref.  10).  The  Pten- null  prostate  cancer  data  were 
generated  in  our  laboratory  and  previously  published  (11).  The  Myc- 
transgenic  prostate  data  were  obtained  online7  (12).  Based  on  UniGene 
mapping,  3,641  unique  genes  were  shared  between  the  MPEDB  and  U74Av2 
microarrays.  Of  these,  3,593  were  in  common  with  the  7,993  genes  used  in 
the  time  course  data.  Significance  analysis  of  microarrays  software  was  used 
to  identify  differentially  expressed  genes  across  prostate  development  or 
between  cancer  and  normal  samples  (13).  Complete  linkage  hierarchical 
clustering  was  performed  using  Cluster  3.0  software  (Eisen  Lab)  without 
weighting. 

Principal  component  analysis  (PCA)  was  performed  using  R.  Expression 
values  for  each  of  the  7,993  genes  were  first  centered  to  mean  zero  across 
the  developmental  time  points  spanning  E15.5-P90.  We  performed  PCA  for 
all  7,993  mean-centered  genes.  The  percent  temporal  variance  captured  by 
each  of  the  first  five  temporal  PCs  was  29.8%,  16.3%,  8.8%,  7.6%,  and  5.0%, 
with  100%  of  the  variance  captured  by  the  first  41  PCs.  To  project  the  mouse 
prostate  cancer  models  onto  the  development  space,  we  implemented  PCA 
on  single-channel  expression  values  for  the  mouse  prostate  development 
data  and  followed  the  methods  described  by  Kho  and  colleagues  (14). 

Gene  Ontology  analysis.  Functional  analysis  of  the  genes  was  carried 
out  using  a  Bayesian  integrated  analysis  (see  Supplementary  Data  for 
detailed  methods).  From  the  Gene  Ontology  (GO)  database,  we  obtained 
data  on  208  biological  processes,  64  cellular  components,  and  151  molecular 
functions.  These  were  selected  from  levels  2  to  6  of  the  GO  hierarchy  that 
were  represented  by  >10  distinct  genes  on  the  array  to  enable  robust 
statistical  analysis. 

Analyses  of  human  prostate  cancer  gene  expression  and  outcomes. 

The  human  orthologs  for  the  91  unique  branching  morphogenesis  Genbank 
identifiers  were  identified  using  the  HomoloGene  database  at  the  National 
Center  for  Biotechnology  Information  and  mapped  to  each  of  three 
independent  prostate  cancer  data  sets  (see  Supplementary  Data).  For  each 
gene,  expression  differences  between  normal  and  cancer  were  assessed  by 
paired,  two-sample  t  test  of  log2-transformed  signal  values.  To  evaluate  the 
power  for  branching  morphogenesis  genes  to  predict  prostate  tumor 
behavior,  we  evaluated  human  prostate  cancer  gene  expression  data  sets 
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with  attendant  clinical  outcomes  (see  Supplementary  Data).  Logistic 
regression  was  then  used  to  measure  the  association  between  prostate- 
specific  antigen  (PSA)  relapse  and  gene  expression  levels.  We  selected  the 
genes  associated  with  recurrent  disease  at  the  P  =  0.05  level  and  used  them 
to  classify  the  tumor  samples  into  two  groups  with  k-means  clustering  (see 
Supplementary  Data  for  methods).  The  difference  in  disease-free  survival, 
defined  as  no  PSA  relapse  >0.1  in  any  follow-up  time  point,  for  the  two 
groups  was  quantified  using  the  survDiff  function  in  R. 


Results 

Patterns  of  gene  expression  associate  with  stages  of  prostate 
development:  prostate  induction,  branching  morphogenesis, 
and  secretory  differentiation.  The  major  events  of  mouse 
prostate  development  can  be  summarized  in  three  steps:  (a) 
prostate  induction  at  E17.5,  ( b )  branching  morphogenesis  at  E18.5 
through  postnatal  day  15,  and  (c)  secretory  differentiation  at 
postnatal  days  25  to  45.  To  characterize  the  transcriptional 
program  associating  with  these  events,  we  profiled  gene  expression 
in  the  prostate  or  prostate  precursor  tissues  at  seven  time  points 
corresponding  to  key  stages  of  prostate  organogenesis:  embryonic 
days  E15.5,  E16.5,  E17.5,  and  E18.5  and  postnatal  days  P7,  P30,  and 
P90  (Fig.  1).  We  generated  three  independent  biological  samples  for 
each  time  point  and  each  sample  was  hybridized  in  a  replicate 
design  to  a  cDNA  microarray  designed  to  assess  gene  expression  in 
the  mouse  prostate  gland  (10,  15).  To  make  the  time  points  directly 
comparable,  a  common  reference  RNA  consisting  of  E14.5  male 
UGS  was  included.  Comparing  transcript  levels  derived  from 
different  developmental  stages  back  to  E14.5  reflects  the  unfolding 
program  of  prostate  development  in  relationship  to  the  most 
undifferentiated  state.  Overall,  70%  of  the  expressed  genes  (5,586 
unique  transcripts)  were  significantly  changed  over  the  time  course 
[false  discovery  rate  (FDR)  <  1%].  Because  of  the  large  number  of 
replicates  (n  =  6  per  time  point),  we  were  able  to  statistically  detect 
small  changes  in  gene  expression  but  chose  to  arbitrarily  confine 
further  analyses  to  the  significant  genes  with  the  highest  temporal 
variance  during  organ  development  (mean  of  5.6-fold  difference; 
range,  2.1-fold  to  237.9-fold;  n  =  2,000). 

Hierarchical  clustering  of  genes  with  the  highest  variance  over 
the  developmental  time  course  revealed  several  distinct  profiles 
(Fig.  1).  One  pattern,  composed  of  371  genes,  was  coincident  with 
prostate  induction,  in  which  expression  peaked  in  E16.5  or  E17.5 
UGS,  decreased  at  E18.5  and  day  7,  and  then  fell  sharply  through 
puberty  and  adulthood.  This  cluster,  designated  the  prostate 
inducer  profile,  included  secreted  frizzled-related  proteins  1  and  2 
(Sfrpl  and  Sfrp2)  and  adult  hemoglobin  a  and  p  chains  ( Hba-al 
and  Hbb-bl ).  Secreted  frizzled-related  proteins  are  a  class  of  Wnt 
pathway  inhibitors  that  modulate  Wnt  signaling  through  direct 
binding  to  Wnt  ligands  as  well  as  binding  to  frizzled  receptors  (16). 
Sfrpl  promotes  prostate  growth  and  branching  morphogenesis  in 
cultures  of  rat  neonatal  ventral  prostate  and  stimulates  the  growth 
of  immortalized  human  prostate  epithelium  (17).  To  our  knowl¬ 
edge,  the  functional  role  of  Sfrp2  in  prostate  induction  has  not  been 
characterized.  Hba-al  and  Hbb-bl  encode  the  major  hemoglobins 
of  adult  erythroid  cells.  Embryonic  forms  of  hemoglobin  (Hba-x, 
Hbb-y,  and  Hbb-z)  were  detectable  at  extremely  low  levels  in  UGS 
before  prostate  budding  and  not  expressed  at  subsequent  time 
points  (data  not  shown).  The  expression  of  hemoglobin  genes  in 
nonerythroid  tissue  was  unexpected  but  not  unprecedented  (18). 
Of  interest,  a  transgenic  mouse  model  designed  to  study 
erythroleukemia  using  fetal  y-globin  or  [Vglobin  promoters  to 
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Figure  1.  Gene  expression  programs 
across  temporal  stages  of  prostate 
development.  Hierarchical  clustering  of  the 
2,000  unique  genes  with  the  highest 
temporal  variance  across  mouse  prostate 
development  is  depicted  as  a  heat  map  in 
which  red  indicates  higher  expression 
relative  to  E14.5  UGS  and  green  indicates 
lower  expression.  E15.5,  E16.5,  E17.5,  and 
El 8.5  are  samples  of  embryonic  UGS, 
whereas  days  7,  30,  and  90  are  of 
postnatal  prostate.  Cluster  dendrograms 
show  that  repeated  samples  from  the  same 
time  point  cluster  together,  excepting  that 
the  closely  related  time  points  of  El  6.5  and 
El 7.5  UGS  are  not  distinguished. 

The  prostate  inducer  profile,  branching 
morphogenesis  profile,  and  secretory 
differentiation  profile  are  identified  to  the 
right  of  the  heat  map  with  a  black,  red, 
and  turquoise  bar,  respectively.  Examples 
of  specific  genes  exhibiting  these  profiles 
are  given  as  HUGO  identifiers  along  with 
plots  of  their  behavior  over  the  time  course. 
The  Y  axes  are  log2  ratios  relative  to 
E14.5  UGS. 
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drive  SV40T  antigen  was  unexpectedly  found  to  develop  prostate 
carcinoma  (19). 

To  validate  the  temporal  profile  of  Sfrpl,  Sfrp2,  Hba-al,  and  Hbb- 
bl  expression,  we  performed  quantitative  RT-PCR  at  each  of  the 
seven  time  points,  producing  results  that  were  highly  concordant 
with  the  microarray  measurements  (Supplementary  Fig.  S2). 
Immunoblots  with  antiserum  recognizing  adult  hemoglobins 
confirmed  high  protein  levels  at  E16.5  with  a  rapid  decline  over 
subsequent  time  points  (Supplementary  Fig.  S3),  and  immunoflu¬ 
orescence  localized  hemoglobin  protein  expression  to  the  UGS 
mesenchyme,  with  no  evidence  of  expression  in  the  epithelial  cells 
of  the  developing  prostate  (Supplementary  Fig.  S4). 

We  expected  genes  associated  with  branching  morphogenesis  to 
increase  expression  levels  after  prostate  induction  (E17.5),  peak  at 
day  7,  and  fall  during  puberty  and  adulthood.  Of  the  most  variant 
genes,  108  (5.4%)  strongly  fit  this  profile,  including  platelet-derived 
growth  factor  a.  {Pdgfa),  follistatin  {Fst),  bone  morphogenetic 
protein  1  {Bmpl),  inhibin  ct  {Inha),  activin  A  receptor  type  II-like 
1  {Acvrl),  and  TGFfl  receptor  3  ( Tgfbr3\  Fig.  1).  Bmpl  and  Inha  are 
secreted  morphogens  in  the  TGF(i  superfamily,  whereas  Fst  is  a 
secreted  inhibitor  of  the  TGF  [J>  family  member  activin  A.  The  high 
representation  of  ligands  (Bmpl,  Inha,  and  Fst)  and  receptors 


(Acvrl  and  Tgfbr3)  that  modulate  TGFfi  signaling  suggests  a  role 
for  the  TGF  (3  pathway  specifically  in  prostate  development.  This  is 
supported  by  evidence  that  TGF[J»,  activin  A,  and  Fst  each  influence 
prostate  branching  morphogenesis  in  neonatal  rodent  ventral 
prostate  cultures  (20,  21). 

The  third  profile  we  examined  were  genes  associated  with 
terminal  secretory  differentiation.  The  mouse  prostate  begins  to 
produce  secretory  proteins  that  contribute  to  the  seminal  fluid  just 
before  puberty,  at  —  30  days  postnatal.  We  expected  these  genes  to 
exhibit  low  or  absent  expression  at  early  time  points,  rise  at  ~  day 
30,  and  remain  high  after  maturity  (22).  Three  hundred  and  forty- 
two  of  the  temporally  variable  genes  in  our  analysis  fit  this  profile 
(Fig.  1).  All  of  the  known  mouse  prostate  secretory  proteins  were 
present  in  this  group,  including  probasin  {Pbsn ),  spermine-binding 
protein  {Sbp),  serine  protease  inhibitor  Kazal  type  3  {Spink3), 
transglutaminase  4  (Tgm4),  seminal  vesicle  secretion  2  {Svs2),  and 
prostate  (3  defensin  1  {Pbdl/Defb37). 

A  previous  report  by  Abbott  and  colleagues  (15)  used  abundance 
measurements  of  expressed  sequence  tags  determined  from 
developmental  stage-specific  cDNA  libraries  to  identify  285  genes 
altered  during  mouse  prostate  genesis.  Of  these  genes,  192  were 
present  on  the  microarrays  used  in  the  present  study.  We  compared 
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the  “virtual”  expression  analysis  of  these  genes  with  transcript 
levels  measured  by  microarray  hybridization  over  roughly  the  same 
time  course.  Hierarchical  clustering  revealed  concordance  between 
the  previous  in  silico  analysis  and  the  microarray  results,  with  most 
genes  exhibiting  qualitatively  similar  temporal  profiles  (data  not 
shown). 

Gene  functions  associated  with  specific  stages  of  mouse 
prostate  development.  To  determine  how  developmental  shifts  in 
mouse  prostate  gene  expression  correlate  with  biological  function, 
we  mapped  GO  identifiers  for  biological  processes,  molecular 
functions,  and  cellular  components  onto  genes  represented  on  the 
mouse  prostate  microarray.  We  used  a  Bayesian  model  that 
incorporates  uncertainty  associated  with  microarray  normalization 
and  gene  classification  to  assess  the  probability  of  functional 
enrichment  (23).  We  confined  our  analysis  to  GO  terms  that  were 
represented  by  at  least  15  genes  to  ensure  that  functional 
conclusions  were  not  drawn  from  categories  with  little  overall 
representation.  We  found  that  several  broad  functional  classes, 
such  as  morphogenesis  (00:0009653),  cell  communication 
(00:0007154),  and  cell  proliferation  (00:0008283),  were  active  at 
earlier  time  points  (P  >  0.9),  whereas  immune  response 
(00:0006955)  and  transporter  activity  (00:0005215)  were  enriched 
in  puberty  and  adulthood. 

To  examine  functional  activities  specifically  associated  with 
prostate  induction,  branching  morphogenesis,  and  secretory 
differentiation,  we  estimated  the  proportion  of  up-regulated  genes 
assigned  a  GO  term  at  each  time  point  and  compared  this  with  the 
proportion  on  the  microarray  as  a  whole.  Hierarchical  clustering  of 
422  GO  categories  revealed  functions  that  closely  followed  the 
prostate  inducer  and  branching  morphogenesis  profile,  but 
surprisingly,  no  functions  strongly  fit  the  secretory  differentiation 
profile  (Supplementary  Fig.  SI).  Although  many  GO  categories 
peaked  at  days  30  or  90,  most  of  these  also  exhibited  a  second 
minor  peak  at  E16.5  or  E17.5  (Supplementary  Fig.  SI).  Functions 
associated  with  the  prostate  inducer  profile  included  cell-cell 
signaling  (00:0007267),  cell-cell  adhesion  (00:0016337),  cell  cycle 
(G0:0007049),  and  the  Wnt  signaling  pathway  (00:0016055; 
Supplementary  Fig.  SI).  In  addition  to  Sfrpl  and  Sfrp2,  we  found 
that  expression  of  Wnt2,  Wnt4,  Wnt5a,  Dkk2,  and  13  other  Wnt 
pathway  components  (out  of  33  genes  on  the  array  annotated  to 
00:0016055)  followed  the  prostate  inducer  profile.  The  cell-matrix 
adhesion  (00:0007160)  and  integrin-mediated  signaling 
(00:0007229)  closely  followed  the  branching  morphogenesis 
profile.  Integrin  signaling  is  crucial  for  ductal  branching  morpho¬ 
genesis  and  several  integrins  are  developmentally  regulated 
(24,  25).  We  found  that  integrins  a3,  a4,  a6,  a7,  (31,  (34,  and  (36 
as  well  as  integrin-linked  kinase  and  RhoA  peaked  in  expression  at 
postnatal  day  7,  during  the  time  of  active  branching  morphogen¬ 
esis.  Taken  together,  these  analyses  suggest  involvement  of  the  Wnt 
pathway  in  prostate  induction,  whereas  integrins  and  TGF(3 
signaling  may  play  a  more  prominent  role  in  branching 
morphogenesis. 

Genetically  engineered  mouse  models  of  prostate  cancer 
share  molecular  features  with  early  prostate  development.  It 

has  been  proposed  that  common  mechanisms  are  shared  between 
development  and  cancer  in  which  neoplastic  lesions  recapitulate 
the  events  of  normal  development  in  reverse,  progressing  from  a 
differentiated  to  a  dedifferentiated  state  (26).  Importantly,  key 
features  of  organogenesis  involve  cellular  characteristics  that  are 
also  hallmarks  of  neoplasia.  In  this  context,  expression  profiling  of 
mouse  cerebellum,  lung,  and  colon  has  revealed  strong  organ- 


specific  connections  between  developmental  programs  and  cancers 
arising  in  these  tissues  at  a  genome-wide  level  (14,  27).  To  explore 
relationships  between  prostate  development  and  neoplasia,  we 
evaluated  two  well-characterized  genetically  engineered  mouse 
models  of  prostate  carcinogenesis,  a  prostate-specific  deletion  of 
the  Pten  tumor  suppressor  gene  (11),  and  a  prostate-specific 
overexpression  of  the  Myc  oncogene  (12).  Mice  with  homozygous 
deletions  of  Pten  develop  invasive  prostate  cancer  by  9  weeks  of 
age.  A  previous  study  using  microarray-based  quantitation  of  gene 
expression  in  Pten _/~  tumors  identified  285  up-regulated  and  241 
down-regulated  genes  (>2-fold  difference;  FDR  <  15%)  compared 
with  normal  prostates  from  litter-matched  controls  (11).  To 
determine  relationships  between  genes  regulated  during  prostate 
oncogenesis  and  the  developmental  phases  of  the  prostate,  we 
examined  the  behavior  of  the  526  Pten^-  tumor-associated  genes 
over  the  prostate  development  time  course  by  PCA  using  an 
approach  described  by  Kho  and  colleagues  (14).  In  normal 
development,  the  dominant  overall  trend  involved  genes  expressed 
highly  at  either  end  of  the  developmental  spectrum.  This  pattern 
was  reflected  by  the  first  temporal  PC  (PCI),  where  those  genes 
expressed  at  high  levels  early  in  development  and  subsequently 
monotonically  decrease  to  low  levels  at  prostate  maturity  receive  a 
negative  first  PC  (PCI  <  0)  and  genes  expressed  at  low  levels  in 
embryogenesis  and  high  levels  at  maturity  comprise  a  positive  first 
PC  (PCI  >  0).  Using  the  convention  of  Kho  and  colleagues,  we 
labeled  genes  with  negative  PCI  values  the  prostate  early  mouse 
partition  (PEMP)  and  the  cohort  with  positive  PCI  values  the 
prostate  late  mouse  partition  (PLMP;  Fig.  2 A).  Of  the  up-regulated 
genes  in  Pten  tumors,  143  (50.1%)  mapped  to  negative  PCI 
coordinates,  partitioning  according  to  early  prostate  development, 
and  142  (49.9%)  mapped  to  positive  PCI  locations.  Among  the  241 
genes  down-regulated  in  Pten tumors,  32  (13.3%)  segregated 
according  to  the  PEMP  and  209  (86.7%)  with  late  development 
PLMP  (Fig.  2 B).  To  address  how  many  genes  we  would  expect  to 
segregate  to  early  or  late  development  by  chance  alone,  we  selected 
526  genes  from  the  developmental  time  course  at  random  and 
performed  PCA  analysis.  Repeated  random  samples  of  526  genes 
showed  that  44.3%  and  55.7%  segregated  to  early  and  late 
development,  respectively.  Therefore,  more  PEMP  genes  are 
expressed  highly  in  Pten  tumors  than  expected  by  chance 
(50.1%  compared  with  44.3%)  and  more  PLMP  genes  are  down- 
regulated  (86.7%  compared  with  55.7%).  Hypothesis  testing  for 
equality  of  proportions  revealed  that  the  enrichment  of  both  the 
up-regulated  Pteri  tumor  genes  in  the  expression  program  of 
early  prostate  development  and  down-regulated  genes  in  late 
development  was  highly  statistically  significant  ( P  <  10~5).  The 
odds  ratio  of  an  up-regulated  Pten  tumor  gene  segregating  with 
early  development  was  6.6  (95%  confidence  interval,  4.2-10.2;  P  < 

itr5,  x2). 

To  investigate  if  the  association  we  observed  between  prostate 
cancer  and  prostate  development  is  generalizable,  we  analyzed  a 
second  model  of  prostate  carcinogenesis  generated  by  over¬ 
expressing  the  Myc  oncogene  in  prostate  epithelium.  Microarray- 
based  expression  profiling  studies  of  Myc-overexpressing  prostate 
tumors  were  previously  reported  (12),  and  we  identified  165 
differentially  expressed  genes  (>2-foId  differential  expression;  FDR 
<  15%)  from  this  study  with  corresponding  features  in  our 
developmental  time  course  experiments.  PCA  analysis  of  the 
behavior  of  these  genes  in  the  developmental  time  course  revealed 
that  74.2%  (49  of  66)  of  Myc  up-regulated  genes  segregated  to  early 
development,  whereas  87.9%  (87  of  99)  of  Myc  down-regulated 
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Figure  2.  Gene  expression  maps 
associating  prostate  cancer  and  early 
development.  A,  an  “egg  plot”  depicts 
the  7,993  genes  expressed  over  the 
developmental  time  course  plotted 
according  to  PCI  and  PC2,  where  PCI 
captures  39%  of  the  variance.  Genes  with  a 
negative  PCI  tend  to  steadily  increase  over 
developmental  time  (early  genes,  gray 
dots),  whereas  genes  with  a  positive  PCI 
tend  to  steadily  increase  over  the  time 
course  (late  genes,  black  dots).  For 
comparison,  heat  maps  depict  the 
developmental  behavior  of  the  early  genes 
(with  the  most  negative  PCI),  and  the 
late  genes  (with  the  most  positive  PCI) 
above  the  corresponding  dots  on  the  plot, 
where  dark  indicates  greater  relative 
expression  and  light  indicates  lower 
expression.  The  seven  time  points  are 
ordered  on  the  heat  maps  from  El  5.5 
on  the  left  through  adult  on  the  right 
(E15.5  UGS,  E16.5  UGS,  E17.5  UGS, 

El 8.5  UGS,  day  7  prostate,  day  30 
prostate,  and  day  90  prostate).  B,  prostate 
developmental  expression  profile  of 
562  genes  that  are  regulated  2-fold 
(FDR  <  15%)  in  Pfen-null  tumors  compared 
with  wild-type  tissue.  Gray  triangles, 

Pten  down-regulated  genes;  black  circles, 
up-regulated  genes.  C,  prostate 
developmental  expression  profile  of 
165  genes  that  are  regulated  2-fold 
(FDR  <  15%)  in  prostates  overexpressing 
Myc  compared  with  wild-type  tissue.  Gray 
triangles,  Myc  down-regulated  genes; 
black  circles,  up-regulated  genes.  Note 
the  correlation  of  prostate  cancer 
up-regulated  genes  (black  circles)  with 
early  developmental  genes  and  cancer 
down-regulated  genes  (gray  triangles) 
with  late  developmental  genes. 
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genes  were  associated  with  late  development  (Fig.  2C).  Proportions 
testing  indicated  that  these  enrichments  were  highly  significant 
(P  <  1(T5).  The  odds  ratio  for  Myc  up-regulated  genes  being 
associated  with  early  prostate  development  was  20.9  (95% 
confidence  interval,  9.2-47.3;  P  <  10  ~5,  x.2)- 

To  determine  if  the  association  between  prostate  cancer  gene 
expression  and  aspects  of  normal  organogenesis  simply  represents 
a  generic  developmental  profile,  rather  than  organ-specific 
developmental  states,  we  analyzed  the  prostate  cancer  expression 


signatures  in  the  context  of  a  temporal  gene  expression  profile  of 
mouse  lung  development,  an  organ  system  that  also  involves 
branching  morphogenesis  (28).  Using  the  PCA  used  for  the  prostate 
studies,  we  found  no  significant  associations  with  stages  of  lung 
morphogenesis. 

Genes  altered  in  murine  prostate  adenocarcinoma  map  to 
the  branching  morphogenesis  stage  of  prostate  development. 

We  next  sought  to  place  the  cellular  phenotype  of  neoplastic 
prostate  epithelium,  represented  by  its  molecular  program  of 
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expressed  genes,  on  the  continuum  of  prostate  development.  We 
first  specifically  evaluated  associations  between  Pten  and  Myc 
cancer  profiles  with  the  induction,  branching  morphogenesis,  and 
secretory  differentiation  profiles  by  comparing  the  proportion  of 
developmental  profile-related  genes  that  received  positive 
t  statistics  (up-regulated)  in  the  tumors  to  the  expected  proportion 
based  on  the  number  of  genes  with  positive  t  statistics  in  the  entire 
data  set.  The  prostate  inducer  profile  was  significantly  enriched  in 
the  Pten  tumors  but  not  in  the  Myc-driven  tumors  [Pten  :  193 
of  294  (65.6%)  genes  up-regulated,  P  <  0.00001,  proportions  test;  Myc: 
113  of  193  (58.5%)  genes  up-regulated,  P  =  0.27,  proportions  test]. 
The  branching  morphogenesis  profile  was  strongly  enriched  in  both 
the  Pten~ /_  and  Myc-driven  cancers  [Pten:  72  of  91  (79.1%)  genes  up- 
regulated,  P  <  0.00001,  proportions  test;  Myc:  49  of  64  (76.6%)  genes 
up-regulated,  P  <  0.0005,  proportions  test].  Conversely,  over  70%  of 
genes  associated  with  the  secretory  differentiation  profile  were 
down-regulated  in  both  Pten- null  and  Myc-driven  cancers.  We  next 
projected  the  expression  profiles  generated  from  individual  Myc  or 
Pten~/ ~  tumors  onto  the  genomic  developmental  trajectory  of  the 
mouse  prostate  represented  by  PCs  (Fig.  3/1 ).  Pten  tumors  localized 
in  a  tight  cluster  between  P7  and  P30,  whereas  Myc  tumors  were 
more  dispersed  over  developmental  space  slightly  preceding  and 
slightly  following  postnatal  day  7  (Fig.  35),  a  finding  that  may  reflect 
the  pleiotropic  activities  attributable  to  the  Myc  protein  (29).  We  are 
not  aware  of  any  morphologic  features  differing  between  these 
models  that  would  account  for  their  placement  on  the  develop¬ 
mental  time  course.  Together,  these  results  are  concordant  with  our 
findings  that  up-regulated  and  down-regulated  genes  in  murine 
prostate  cancer  segregate  with  early  and  late  development, 
respectively.  In  addition,  these  data  suggest  that  genes  comprising 
the  branching  morphogenesis  profile  may  contribute  to  processes 
influencing  carcinogenesis. 

Genes  comprising  the  branching  morphogenesis  program 
are  altered  in  human  prostate  carcinoma.  The  complex 


developmental  process  of  branching  morphogenesis  involves 
several  features  that  are  also  operative  in  invasive  prostate  cancers. 
These  include  cellular  processes  contributing  to  cell  movement, 
adhesion,  invasion,  division,  and  death  that  are  extensively 
influenced  through  interactions  with  surrounding  mesenchyme 
and  extracellular  matrix.  To  determine  if  genes  operative  in  the 
normal  branching  morphogenesis  program  exhibit  elevated  activity 
in  human  prostate  cancers,  we  evaluated  transcript  abundance 
levels  of  orthologous  human  genes  measured  by  microarray-based 
profiling  studies  of  human  tumors.  These  comprised  two  studies, 
True  and  colleagues  (30)  and  Tomlins  and  colleagues  (31),  where 
microdissected  epithelium  was  acquired,  and  one  study  reported 
by  Lapointe  and  colleagues  (32)  that  used  macrodissected  tissue 
samples.  Overall,  in  each  data  set,  ~  20%  of  the  genes  comprising 
the  murine  branching  morphogenesis  program  were  found  to  be 
statistically  increased  in  human  prostate  cancers  relative  to  benign 
tissue.  Although  the  gene  representation  varied  across  studies,  in  a 
large  part  due  to  different  probes  present  on  the  different 
microarray  platforms,  several  genes  exhibited  significant  increases 
in  prostate  cancers  consistently  including  peroxiredoxin  4  ( PRDX4 ), 
SLC43A1/POV1 ,  and  the  DNA  (cytosine-5-)-methyltransferase  3a 
(DNMT3A:  Fig.  4).  Importantly,  many  of  the  genes  comprising  the 
branching  morphogenesis  network  have  not  been  extensively 
studied  in  the  context  of  prostate  neoplasia. 

We  next  sought  to  determine  if  genes  associated  with  branching 
morphogenesis  were  associated  with  malignant  tumor  character¬ 
istics.  We  evaluated  a  data  set  reported  by  Stephenson  and 
colleagues  (33)  that  used  microarrays  to  profile  transcript  levels  in 
prostate  tumors  from  79  patients  with  attendant  clinical  follow-up 
delineating  biochemical  relapse  after  radical  prostatectomy.  Of  the 
91  genes  comprising  the  BMP,  84  had  orthologs  represented  in  this 
data  set.  A  single  estimate  of  the  expression  level  in  each  sample 
was  generated  and  the  association  between  PSA  relapse  and 
transcript  levels  was  determined  by  logistic  regression.  We  found 
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Figure  4.  Genes  comprising  the  branching  morphogenesis  signature  are  dysregulated  in  human  prostate  cancer  and  associate  with  outcome.  Microarray-based 
measurements  of  SLC43A1  (A),  PRDX4  (B),  and  DNMT3A  (C)  transcript  abundance  levels  between  normal  (A/;  circles)  and  cancerous  (C)  epithelium  ( triangles ) 
are  plotted  from  three  independent  data  sets.  For  True  and  colleagues  data  (n  =  35),  differences  between  normal  and  cancer  were  assessed  by  paired,  two-sample 
t  test  of  log2-transformed  signal  values.  For  Lapointe  and  colleagues  data  (n  =  41),  differences  between  normal  and  cancer  were  assessed  by  paired,  two-sample 
f  test  of  log2-transformed  ratios  to  cell  line  control.  For  Tomlins  and  colleagues  data  (n  =  31),  differences  between  normal  and  cancer  were  assessed  by  unpaired, 
two-sample  t  test  of  log2-transformed  ratios  to  benign  prostate  pool.  D,  a  14-gene  signature  associated  with  branching  morphogenesis  was  used  to  classify  prostate 
cancer  samples  by  disease  recurrence  as  determined  by  PSA  relapse.  The  difference  in  disease-free  survival,  defined  as  no  PSA  relapse  >0.1  in  any  follow-up 
time  point,  for  the  two  groups  was  quantified  using  the  survDiff  function  in  R  and  the  P  value  was  reported  by  log-rank  test. 


that  14  genes  were  significantly  associated  with  disease  relapse 
(P  <  0.05)  and  a  classifier  using  this  gene  signature  provided 
discriminatory  outcome  information  associating  with  disease 
recurrence  (x2  =  14.4;  P  <  0.001;  Fig.  4).  Further  confirmation  of 


the  independent  predictive  power  of  genes  comprising  the 
branching  morphogenesis  profile  will  require  a  directed  study  with 
long-term  clinical  outcomes  that  incorporates  other  risk  factors 
associated  with  tumor  behavior. 
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Discussion 

Carcinomas  arising  in  a  diverse  range  of  organs  and  cell  types 
are  known  to  display  immature  features  and  are  accompanied  by 
marked  changes  in  gene  expression  (34).  A  subset  of  genes 
aberrantly  expressed  in  tumor  cells  is  known  to  normally  exhibit 
highly  compartmentalized  spatial  and  temporal  expression  pat¬ 
terns  localized  to  specific  stages  of  embryonic  development  (9,  35). 
Several  embryonic  proteins,  such  as  carcinoembryonic  antigen  and 
a-fetoprotein,  have  been  developed  into  useful  diagnostic  tumor 
markers.  However,  with  few  exceptions,  the  overall  relationships 
between  tumors  and  developmental  programs  defined  at  the 
genetic  level  have  not  been  evaluated  (14,  27).  In  this  study,  we 
sought  to  determine  molecular  features  underlying  events  linked  to 
normal  developmental  processes  in  the  prostate  gland  and  those 
accompanying  neoplastic  transformation.  Overall,  we  identified 
significant  correlations  between  gene  expression  profiles  represent¬ 
ing  early  stages  of  prostate  organogenesis  and  two  distinct  mouse 
models  of  prostate  carcinoma  arising  in  the  context  of  Pten  loss  or 
Myc  overexpression.  In  general,  transcripts  differentially  up- 
regulated  in  carcinomas  were  more  likely  to  also  be  expressed 
early  in  development  and  decline  with  progressive  organ 
maturation,  whereas  genes  expressed  highly  in  the  mature 
differentiated  prostate  gland  were  down-regulated  in  carcinomas. 
Similar  overall  trends  were  observed  in  studies  comparing  gut 
development  and  colorectal  tumors  (27)  and  of  cerebellar 
development  and  medulloblastomas  (14).  Together,  these  studies 
provide  strong  support  for  a  general  reactivation  of  primitive 
cellular  programs  operating  to  govern  a  range  of  phenotypic 
opportunities  involving  cell  position,  division,  motility,  and 
invasion  that  are  provided  by  intrinsic  and  extrinsic  cues. 

The  molecular  profiles  derived  from  prostate  cancers  exhibited 
the  greatest  association  with  the  branching  morphogenesis  stage  of 
prostate  development.  The  process  of  branching  morphogenesis 
involves  a  complex  interplay  of  cellular  events  that  in  many  ways 
recapitulates  features  of  malignant  cells  (36).  The  process  is  critical 
for  the  formation  of  arborized  organs  that  span  the  development  of 
tracheal  networks  in  insects  to  a  diverse  array  of  human  tissues 
that  include  the  pancreas,  lung,  salivary  gland,  kidney,  breast,  and 
prostate  (6,  37).  Branching  morphogenesis  entails  reorganization  of 
epithelial  tissues  to  form  complex  but  highly  structured  tubular 
assemblies  that  function  to  produce  and  transport  fluids  and  gases 
over  large  surface  areas  (36).  A  series  of  sequential  and  often 
iterative  biochemical  and  biomechanical  steps  are  required  for  the 
proper  construction  of  the  networks,  and  many  of  the  key 
molecular  features  governing  these  processes  have  been  elucidated. 
A  critical  component  of  the  process  involves  the  invagination  and 
subsequent  invasion  of  epithelial  buds  and  outgrowths  into 
surrounding  mesenchyme.  Importantly,  the  incursion  of  epithelium 
is  highly  dependent  on  signals  derived  from  the  surrounding 
stroma,  an  attribute  increasingly  recognized  to  play  an  important 
role  in  carcinogenesis.  Key  regulators  of  these  processes  include 
FGF10,  Shh,  Bmp4,  TGF(3,  as  well  as  members  of  the  matrix 
metalloproteinase,  a  disintegrin  and  metalloproteinase,  and  serine 
protease  families  (38,  39).  The  matrix  metalloproteinases  are  of 
particular  interest  due  to  their  integral  roles  in  regulating 
epithelial-mesenchymal  cross-talk  and  their  influence  on  migratory 
processes  through  proteolysis  of  matrix  molecules  and  the 
generation  of  motogens  such  as  laminin-5  fragments  (40).  The 
leading  edge  of  migrating  epithelial  cells  has  been  shown  to  exhibit 
mesenchymal  phenotypes,  allowing  for  penetration  through  matrix 


and  stromal  cell  tissue  constituents,  a  process  also  well  described 
in  the  context  of  tumor  cell  epithelium  to  mesenchymal  transition 
(41,  42).  Thus,  tumors  exhibiting  transcript  profiles  congruent  with 
a  branching  morphogenesis  developmental  stage  might  be 
expected  to  exhibit  characteristics  of  enhanced  invasion,  metasta¬ 
sis,  and  early  relapse  after  surgical  resection. 

The  comparative  profiling  studies  we  report  here  identified 
several  genes  that  have  not  been  extensively  studied  in  the 
context  of  normal  developmental  processes  or  prostate  carcino¬ 
ma.  PRDX4  is  a  member  of  a  multifunctional  antioxidant  protein 
family  that  primarily  serves  to  provide  cellular  protection  against 
oxidative  stress  (43).  Peroxiredoxin  family  members  also  regulate 
proliferation,  in  part  through  intracellular  signaling  cascades  that 
apply  hydrogen  peroxide  as  a  molecular  second  messenger.  No 
known  roles  for  PRDX4  have  been  described  in  the  context  of 
development,  although  other  peroxiredoxin  family  members 
interact  with  the  androgen  receptor  (AR)  and  modulate  AR- 
mediated  signaling  (44).  All  of  the  prostate  cancer  studies  we 
evaluated  showed  elevated  expression  of  PRDX4  in  neoplastic 
lesions.  DNMT3A  encodes  an  enzyme  involved  in  de  novo 
methylation  of  genomic  DNA,  a  critical  step  for  regulating 
genomic  imprinting  and  X-chromosome  inactivation.  Prior  studies 
have  shown  aberrant  de  novo  methylation  of  growth-regulatory 
genes  in  human  tumorigenesis  (45).  To  date,  no  specific  role  for 
DNMT3A  has  been  shown  for  genitourinary  tract  developmental 
processes,  although  a  myriad  of  developmental  abnormalities 
results  from  deleting  Dnmt3a  in  mouse  models  and  conditional 
Dnmt3a  mutant  males  show  impaired  spermatogenesis  (46). 
Studies  in  prostate  cancer  have  shown  increased  Dnmt3a 
expression  in  tumors  developing  in  the  TRAMP  model  (47),  and 
increased  DNMT3A  expression  has  been  associated  with  progres¬ 
sion  to  androgen-independent  growth  in  vitro  (47,  48).  SLC43A1 
encodes  a  protein  that  functions  in  the  sodium-independent 
transport  of  neutral  amino  acids.  SLC43A1  was  originally 
identified  as  a  transcript  up-regulated  in  a  clinically  aggressive 
prostate  cancer  and  designated  POV1  (49).  To  date,  there  are  no 
studies  showing  a  mechanistic  role  for  SLC43A1  in  developmental 
processes. 

In  summary,  global  assessments  of  gene  use  in  prostate  cancers 
and  normal  development  provide  intriguing  links  between  the  two 
processes.  This  perspective  allows  for  the  identification  of  specific 
genes  as  well  as  regulatory  patterns,  pathways,  and  networks  that 
operate  to  direct  the  complex  processes  required  for  both  the 
homeostasis  and  evolution  of  normal  and  malignant  tissues. 
Exploiting  normal  developmental  systems  may  provide  a  conve¬ 
nient  and  tractable  model  to  study  mechanisms  that  play 
influential  roles  in  invasive  neoplastic  growth. 
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